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In this paper we study Inozemtsev's su(r7i) quantum spin model with hyperbolic interactions and 
the associated spin chain of Haldane-Shastry type introduced by Frahm and Inozemtsev. We com- 
pute the spectrum of Inozemtsev's model, and use this result and the freezing trick to derive a simple 
analytic expression for the partition function of the Frahm-Inozemtsev chain. We show that the 
energy levels of the latter chain can be written in terms of the usual motifs for the Haldane-Shastry 
chain, although with a different dispersion relation. The formula for the partition function is used 
to analyze the behavior of the level density and the distribution of spacings between consecutive un- 
folded levels. We discuss the relevance of our results in connection with two well-known conjectures 
in quantum chaos. 



I. INTRODUCTION 



o 

^j Over the last few years a significant amount of effort has been devoted to the study of spin chains 

^^ of Haldane-Shastry (HS) type, due to their remarkable integrability properties and their interest in 

^^ connection with several important conjectures in quantum chaos. This class of chains, intimately 

related to integrable dynamical models of Calogero-Sutherland type [IHS], are characterized by 
the fact that the interactions between the spins are both long-ranged and position-dependent. 
O , For instance, in the original Haldane-Shastry chain [H [S] the spins occupy equidistant positions 

I on a circle, the strength of the interactions being inversely proportional to the square of the 

'^ distance between the spins measured along the chord. Historically, this model was introduced 

Cd while searching for a spin chain whose exact ground state coincided with Gutzwiller's variational 

C wavefunction for the one-dimensional Hubbard model in the limit of large on-site interaction [B]- 

I— I [5]. In fact, the HS chain can be obtained in this limit from the Hubbard model with long- 

range hopping in the half-filling regime ^. In Haldane's original paper [4], the spectrum of 
the HS chain with spin 1/2 was inferred on the basis of numerical calculations. In particular, 
it was observed that the levels are highly degenerate, which suggests the presence of a large 
QQ underlying symmetry group and the possible integrability of the model. This symmetry group was 

■^ subsequently identified [10] as the Yangian 3^(slm), where m is the number of internal degrees of 

f^ freedom. As to the model's integrability, it was established around the same time by Fowler and 

ly-^ Minahan |llj using Polychronakos's exchange operator formalism [12^. 

f*""^ The rigorous derivation of the spectrum of the HS chain with arbitrary su(to) spin was carried out 

f^ by Bernard et al. [ISj by taking advantage of its connection with the generalization of Sutherland's 

^-H model to particles with spin [14j . At the heart of this connection is the mechanism known as 

^ Polychronakos's "freezing trick" [TS]. The physical idea behind this mechanism is that when the 

coupling constant of the spin Sutherland (trigonometric) model tends to infinity, the particles 

concentrate around the equilibrium of the scalar part of the potential, so that the dynamical and 

internal degrees of freedom decouple. It can be shown that the coordinates of this equilibrium are 

essentially the HS chain sites, and that in this limit the internal degrees of freedom are governed 

by the chain's Hamiltonian. The freezing trick can also be applied to the spin Calogero (rational) 

model [16], obtaining in this way a spin chain — the so-called Polychronakos-Frahm (PF) chain fT5| 

117) — with non-equidistant sites given by the zeros of the A^-th degree Hcrmite polynomial, N being 

the number of spins. 

The original Calogero and Sutherland models mentioned above are both based on the root 
system yl^v-i, in the sense that the interaction between the particles depends only on their relative 
distance. As shown by Olshanetsky and Perelomov [TH], there are integrable generalizations of these 
models associated with any (extended) root system, the rank of the root system basically coinciding 
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with the number of particles. For this reason, the most studied models of Calogcro-Sutherland (CS) 
type are by far those associated with the BCn root system (including the hyperbolic Sutherland 
model) [THH22] , and to a lesser extent with the Dm system [53J |21] ■ A hyperbolic variant of the 
Sutherland model of A^-i type with an external confining potential of Morse type has also been 
considered in the literature, in both the scalar |23 and the spin [23] cases. 

For all the spin CS models mentioned in the previous paragraphs, a corresponding spin chain of 
Haldane-Shastry type has been constructed by means of the freezing trick [19l HH [23l [24l |27l [28] . 
In the case of the rational and trigonometric chains, this mechanism has been applied to derive a 
closed-form expression for the partition function in terms of the quotient of the partition functions 
of the corresponding spin and scalar dynamical models [22H2^ E51 - IST] , whose spectrum can be 
easily computed. Expanding the partition function in powers of g = e"^'''"'^"^^ one can compute 
the chain's spectrum for relatively large values of N and determine some of its statistical properties. 
A common feature of all of these chains is the fact that when the number of sites is sufficiently 
large the level density is approximately Gaussian. This result, for which there is ample numerical 
evidence, has also been rigorously established in some cases f32]. The knowledge of a continuous 
approximation to the (cumulative) level density is of great importance in the context of quantum 
chaos, as it is used to transform the raw energies so that the resulting "unfolded" spectrum has 
an approximately uniform level density |25j. The distribution of spacings between consecutive 
unfolded levels is widely used for testing the integrable vs. chaotic character of a quantum system. 
Indeed, according to a long-standing conjecture of Berry and Tabor '341 , the spacings distribution 
p{s) of a "generic" quantum system whose classical counterpart is integrable should follow Poisson's 
law p{s) ~ e^^ . On the other hand, the Bohigas-Giannoni-Schmidt conjecture ^Q asserts that 
the spacings distribution of a fully chaotic quantum system is given by Wigner's surmise p{s) = 
{■Ks/2) exp(— 7rs^/4), characteristic of the Gaussian orthogonal ensemble (GOE) in random matrix 
theory. Both of these conjectures have been shown to hold in many different systems, both in 
the integrable ^6l [37] and fully chaotic cases jS^ • Rather surprisingly, the spacings distribution 
of all the integrable spin chains of HS type studied so far is neither of Poisson's nor Wigner's 
type [231 1301 1311 [5HHiD] . Thus, at least in this respect, spin chains of HS type appear to be 
exceptional among the class of integrable systems. 

Unlike their rational or trigonometric counterparts, HS chains with hyperbolic interactions have 
received comparatively less attention in the literature. In particular, the spectrum of the hyperbolic 
BCn chain is not known, whereas that of the Ajv-i chain has been conjectured (27j only for spin 
1/2. The main purpose of this paper is precisely to fill this gap for the hyperbolic chain of A^-i 
type, which we shall henceforth refer to as the Frahm-Inozemtsev (FI) chain. Although at the 
formal level this chain is the hyperbolic analog of the original Haldane-Shastry chain, in practice 
both chains turn out to be quite different. Indeed, while the sites of the HS chain are equidistant on 
a circle, those of the FI chain are not and, moreover, depend on a parameter /3 in a nontrivial way. 
We shall see that the energies of the FI chain also depend on this parameter, which complicates the 
statistical analysis of the spectrum. An important consequence of this dependence is the fact that 
for generic values of /3 the spacings distribution is qualitatively different from that of the integrable 
spin chains of HS type studied so far. At any rate, for all values of (3 the spacings distribution 
of the FI chain is neither of Poisson's nor Wigner's type, in spite of the fact that this chain is 
probably integrable |27j. 

The paper is organized as follows. In Section[ll]we recall the definition of Inozemtsev's hyperbolic 
spin dynamical model and its scalar counterpart, and explicitly construct the Frahm-Inozemtsev 
spin chain by applying the freezing trick to these models. Section [TH] is devoted to the computation 
of the spectrum of Inozemtsev's spin dynamical model, which was partially known only in the case 
of spin 1/2. Our approach is based on relating the Hamiltonian to an auxiliary differential-difference 
operator, which we triangularize by expressing it in terms of suitable Dunkl-Cherednik operators 



of type A [m [55] . In Section IV we use the freezing trick and the results of the previous section 
to derive a closed-form expression for the partition function of the su(to) FI chain. Using this 
expression and some general results for other spin chains of A^-i type [131 SI], we obtain a simple 
formula for the spectrum in terms of the usual motifs jlOj . In particular, this provides a rigorous 
proof of Frahm and Inozemtsev's conjecture for the spectrum in the spin 1/2 case. With the 
help of the partition function, in Section [V] we analyze several statistical properties of the chain's 
spectrum. When /3 = 0{N), our numerical computations show that the level density is Gaussian 
as the number of sites tends to infinity. Taking as the spectrum unfolding function the cumulative 
Gaussian distribution, we have also studied the density of spacings for large N and different values 



of the spin when the parameter /3 is 0{N). Our calculations clearly indicate that the density 
of spacings exhibits the behavior previously found in other chains of HS type only when /3 is an 
integer or a rational with a "small" denominator. The paper ends with two technical Appendices 
in which we prove the existence of a unique solution of the system defining the chain sites, and 
derive a closed-form expression for the mean and variance of the chain's energies. 



II. THE MODELS 

In this section we describe Inozemtsev's hyperbolic spin dynamical model |26| and its corre- 
sponding spin chain of Haldane-Shastry type [27], whose study is the purpose of this paper. The 
Hamiltonian of the spin dynamical model is given by 



ir^J 



sinh [xi ~ Xj ) 



where the sums run from 1 to A^ (as always hereafter, unless otherwise stated), a > 1/2, 6 > 0, 
e = ±1, and the operators Stj permute the spins of the i-th and j-th particles. More precisely, let 
E = (C™)^^ be the space of internal degrees of freedom, and denote by |s) = |si, . . . , sjv), with 
Si € {1, . . . , to}, an element of the canonical basis of S. The action of Sij on |s) is then given by 

Oij\ ■ . ■ , Si, . . . , Sj, . . . ) — \ . . . , Sj , . . . , Si, . . . ) . 

It is well known that the operators Sij can be expressed in terms of the generators f^ of the 
fundamental representation of su{m) for the fc-th particle as 

^ij ~ 2 y, ti tj H J 

^ — ' ■* m 

a=l 

where we have used the normalization tr{t'^tl) = ^6"^. It was shown in Ref. [2S1 that the above 
model is completely integrable for arbitrary m. Regarding its spectrum, in the latter reference only 
a few eigenstates (including the ground state) and their corresponding energies were computed in 
the special case of spin 1/2 (to = 2). The Hamiltonian (fTl) is the spin version of the scalar model 

iJ^'^ = -A + 62 ^(e-2^' - 1)2 + aia - l)J2 sinh"^(a;^ - x^) , (2) 

i i^j 

previously studied in Ref. |25| . In particular, in the latter reference the spectrum of H^'^ was 
computed in closed form, together with the eigenfunctions of the ground state and several excited 
states. It should be noted that, due to the impenetrable nature of the singularities of H and H^'^ in 
the hyperplanes Xi = Xj (z ^ j), the configuration space of both Hamiltonians (II]) and (pi) should 
be taken as one of the Weyl chambers of ^at-i type, say 

C = {x = (xi, . . . , xn) eR^ \xi < X2 < ■■■ < xn} ■ (3) 

We shall next apply Polychronakos's "freezing trick" [T5] to the spin dynamical model (IT]) in 
order to construct its associated spin chain of Haldane-Shastry type. To this end, we rescale the 
strength of the Morse potential as 



and consider the strong coupling limit a 
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where the scalar potential U is given by 

t/(x) = /32 ) '(e-2-' -!)' + >' sinh-"(a;, - x^) . (4) 



We shall show in Appendix |X] that U possesses a minimum in the set C if and only if 

/?>2(7V-1), (5) 

and that this minimum is in fact unique. It follows that, as a -> oo, the eigenfunctions of H 
become sharply peaked around the unique minimum $, = {^i^ . . . , S,n) oi U in C . Since 

H = W+Aah{y.), (6) 

where 

h( "i — - \^ 1 — ebjj , , 

4 ^ sinh^(a;i - Xj) ' 

in the limit a — >■ c» the dynamical and internal degrees of freedom decouple, the latter being 
governed by the Hamiltonian 

^ " MO = ^ E . l^f'% , ■ (8) 

4 ^ smh {^i - Q 



As shown in Appendix [Al the chain sites £,i can be expressed in terms of the zeros (i < ■ ■ ■ < (n 

-0- -- 

2e, = log/3-logCw-,+i. (9) 



of the generalized Laguerre polynomial i^ as 



Recall, in this respect, that Eq. (Is]) is precisely the condition which guarantees that the zeros of 
Lj^ are positive and distinct |45| . In fact, since 14. is invariant under ^j i— > — ^^ + c, with c a 

constant, we can alternatively define the sites of the chain ([8]) by the simpler formula 

C^^^logC- (10) 

We can also express the Hamiltonian (Is]) directly in terms of the zeros Ci of L^~ as 

^ = E77%^(1-^^^.)- (11) 



We shall take Eqs. (p|)-(10l as the definition of the Hamiltonian of the Frahm-Inozemtsev chain, 
although it differs from the normalization used in Ref. [U] by a factor of 2. If e takes the value 1 
(respectively —1) the corresponding FI chain is of ferromagnetic (respectively antiferromagnetic) 
type. We shall sometimes use the more precise notation 'H^ (respectively Ti.^^) to denote the 
Hamiltonian of the ferromagnetic (respectively antiferromagnetic) chain. Note finally that, unlike 
other spin chains of HS type associated with the ^at-i root system, the Hamiltonian of the FI 
chain depends on an essential parameter /3 > 2 (TV — 1) through the zeros C,i. 

III. SPECTRUM OF THE DYNAMICAL MODELS 

In this section we shall compute the point spectrum of Inozemtsev's spin dynamical model M 
and of its scalar counterpart ([2]). As is customary when studying quantum Calogero-Sutherland 
models with spin, we introduce the auxiliary scalar operator 

g = -A + b^ E(e"'^' - 1)^ + ^E ■ ,V '^" . ' (12) 

i I^smh [x^-Xj) 

where Kij is the coordinate permutation operator defined by 

(-l^ijjjy^lj • ■ • J ^ij ■ ■ ■ ,Xj, . . . , Xi^) = j[Xi, . . . , Xj, . . . , Xi, . . . , XjS[) ■ (13) 

The operator H is naturally defined on (a suitable dense subset of) the Hilbert space L^(IR^). On 
the other hand, due to the nature of their singularities on the hyperplanes Xi — Xj (with i ^ j), the 



Hamiltonians H and H^^ act on (appropriate dense subsets of) the Hilbert spaces L'^{C) (g) E and 
L'^{C), respectively. However, proceeding as in Ref. [21], one can show that these Hamihonians 
are isospectral to their extensions H and H^'^ to A(L^(M^) (g) S) and A'^'^(L^(IR^)), respectively, 
where A and A*"^ denote the projectors onto spin and scalar states with any fixed symmetry under 
particle permutations. We shall choose these extensions so that 

This is equivalent to the requirement that 

K,,A = eS,,A , i^,, A^^ = A^^ , 

so that A must project onto spin states with parity e under particle permutations, while A***^ is the 
symmetrizer under coordinate permutations. 

A. Dunkl operators 

In view of the above remarks, in order to compute the point spectra of H and H^'^ it is enough 
to solve the analogous problem for the auxiliary operator H in L'^(M.^). To this end, we introduce 
the gauged Dunkl operators 

2Jr = e"2^' [9^, - 6e^2x, _ ^^ (1 + coth(a;, - Xj))K^j + 6 - ij , (15a) 

2J° = d^^ - 6e"^"^* - a ^ (1 + coth(a;i - Xj))K^j + 2a ^ Kij , (15b) 

in terms of which 

H = -Aj2{j,')'-AbJ2jr+Nb\ (16) 



The operators ( 15 ) are related to the type A Dunkl-Cherednik operators [HI |42] 



Jf =d,^+aJ2 -^- (1 - K^j) , (17a) 



3+ 



■3 



J,° = z,d,^+aY, -^ (1 - K,^) - aJ2^^3 + «(^ - 1) + -^^ (17b) 

by the gauge transformation 

Ji^ ^ PJ^ P^^ : Jz° = pJi P~^ , 

where Zi = e^^' and p is the ground state of the scalar Hamiltonian H^"^, given by 

p = exp[(l + a(iV-l)-fe)^a;,-^^e-2-.] [] | sinh(a;, - x,)!" ■ (18) 

i i i<j 

It is well-known [5(1 that the operators ([It]) preserve the polynomial modules 



T^n = {Y[z'^^ |n, eNo, n, sCn} 



for arbitrary n S No = N U {0}. Proceeding as in Ref. [5S], we seek the eigenf unctions of H in 
pTZn, where n must be chosen so that pTZn C L^(E^). Taking into account that 

nsinh(a;, - x,) oc ^6^'+-^ • [] (e-^-. _ e'^^^) = Y[e^i^''^-^ ■ J] (e"""^- - e-^^') 



we have 



p[| zf • ex [| exp [ (2n, + 2a{N -l) + l-b)x,~ ?^"'"i ' IT N"'"' " ^"'"'T ' 



2 



and thus p Hi -2"' i^ square- integrable provided that 



n,<-{h-l)-a{N-l), l^i^N. 

Hence the maximum value of n such that pTZn C -L^(M^) is given by 

n = max i n' e No | n' < - (6 - 1) - a(iV - 1) [ • (19) 

Note, in particular, that H possesses bound states if and only if the strength of the Morse potential 
satisfies the condition 

b > 2a{N - 1) + 1 . (20) 

B. Triangularization of H 

In order to compute the point spectrum of the auxiliary operator H , we shall construct a non- 
orthonormal basis of pTZn (with n given by Eq. ( |19| ) on which this operator acts triangularly. To 
this end, it suffices to construct a basis of 7?.„ on which the gauge transformed operator 

if'^p-iifp = -4^(jP)'-46^Jr+A^62 (21) 



(cf. Eq. (16)) is represented by a triangular matrix. In fact, the elements of this basis are simply 



the monomials 

</'n = n^r> n=(ni,...,n^)eN^, 

i 

ordered as we shall now explain. Given a multiindex n e N^, we define the associated non- 
increasing multiindex 

[n] = (n,;^ , . . . , ni„ ) , with n^i ^ • • • ^ n^^ . 

For p, p' G [N^], we shall write p -< p' if the first nonzero component of p' — p is positive. In the 
case of two arbitrary multiindices n, n' e N^, we define 

n ^ n' ^==^ [n] -< [n'] . 



Thus, for instance. 



since clearly 



(1,2,3,2)^(1,2,2,6)^(1,1,6,3) 



(3,2,2,1)^(6,2,2,1)^(6,3,1,1) 



Finally, we shall set 0n -< '/'n' if and only if n -< n'. This defines a partial order -< in the set 
of all monomials (/)n S Ti-n- We shall see that the auxiliary operator H' is represented by an 
upper triangular matrix in the basis B of 7?.„ consisting of the monomials </>„, ordered in any way 
consistent with the relation ■<. The proof of this fact can be divided into four steps. 

First of all, it is easy to show that the operators J^^ are strictly upper triangular with respect 
to the basis B, i.e., 

J- </-« - ^ c-„,>„, . (22) 

n'-<n 



Indeed, 






Since the terms in the sum over j ^ i vanish when n^ = rij, the latter sum can be written as 

y. 1 (Z,-/Z.)"'-"^ - 1 y^ 1 (Z,/Z,)"^-"' - 1 

^-^ Zi (Zj/Zi) — 1 ^-^ Zj iZil Zj) — 1 



rii — nj — 1 nj—ni — 1 

E E ^-r'^-^- E E 

j;nj<ni k—0 j',nj>ni k—0 






Thus 



ni— rij — 1 nj—ni — 1 

Ji (f'n — ni4>n-e, + a 2_^ E '?^n-(fe+l)e. + fee3 ^ « 2^ 2^ 0n+fee.-(fc + l)ej j 

j:nj<ni k—0 j\nj>ni k—0 

where e^ is the k-th element of the canonical basis of K^. It can be easily checked that the vectors 
n — ej, n — (fc + l)ei + kej and n + fce^ — (fc + l)ej appearing in the previous expression all precede 
n, which establishes our claim. 

Consider now a non-decreasing niultiindex p G [N^] , and set 

e{pi) = min {j I pj =Pt} , #(pO = card {j \ pj = p^} . 

The second step in the proof of the triangular character of H' with respect to the basis B consists 
in showing that 

J.Vp - A, (p) Cj)p+J2 C^,n'p -^n' , (23) 

n'-;p 



where 



A.(p) =p^ + ^il-b) + a{N + ^ + l- #ip,) - 2£(k)) . (24) 



Indeed, proceeding as before we obtain: 

%pi[j,0-a(iV-l) + i(5-l)-p. 



^ (^.Aj)-i ^^. (^.M)-i ^U^ 

E E (I) - E E (J) -card{,>.|p,=M 

j>i\Pj<Pi k—1 j<i',Pj>Pi k—0 v J/ 

Pi-pj-1 / . fe Pj-Pi-1 / X *: 

Zt \ V ^ X ^ / Zi 



. Zj , 

j>i;Pj<Pi fe=l ^ ' j<i;Pj>Pi fe=l 



E E (f ) - E E 

j<i;Pj>Pi fe=l 

card{j < z I pj > Pj} - card{j > i | p^ = pj . 



Taking into account that 

cardjj <i\pj > pi} = £{pi) - 1, card{j > i \ pj = pj = #(pi) + e{pi) - j - 1 

we have 

Pi-Pj-i Pj-Pi-i 

•^i°'^P = Ai(p)(^p+a ^ ^ (t)p-ke, + ke, - a ^ ^ 0p+fcei-feej , 

j>i\Pj<Pi fc=l 3<i\Pi>Pi fc=l 



which obviously proves our claim. 

Consider next the action of the operator J^ on a basis function ^n with arbitrary n e Nq . A 
computation totally analogous to the previous one shows that 

- j" — n^ — a{N — 1) + - (6 — 1) + a card{j | rij > rii} + a cardjj > i \ rij — ni}\ 0n 

rii — nj — l nj—rii — l 

j<i-^nj<ni j>i;nj>ni j;nj<ni k—1 j;nj>ni fc — 1 

The first two sums in the last expression involve basis functions with multiindices n' such that 
[n'] = [n] , while all the multiindices appearing in the last two sums precede n. Hence we have 

jUn= E ~<n'n<t>n'. (25) 

n';[n']r<[n] 

Although the last identity indicates that the operators J^ need not be triangular with respect 
to the basis B, we shall next show that the sum ^,(J.f)^ appearing in H' is upper triangular in 
the latter basis. This fact, together with Eq. p2| ), implies that H' is represented by an upper 
triangular matrix in the basis B. 

Indeed, given n e N^ let P be any permutation such that n = ^[n], and (with a slight abuse 
of notation) denote also by P the linear operator defined by P(j)n' = (pPn', for all n' € Nq . Since 
{P(f>n'){'2') = (j)n'{P~^z), the operator J^ii'^i*)'^ obviously commutes with P, and therefore 

$:(jp)v„=p5:(j.")>„]. 



Calling [n] = p and using Eqs. (23 1 and ( p5| ) we easily obtain 



E(^^°)'^=E'^°(^^(P)'^P+ E ^n'p'^n' 



J2 A<(p)Vp + J2J2 A,(p)c°„,p0„- + E E c°„,p£?,„.„,0„" , 

i i n'^p i n'^p 

fn"l-<fn'l 



and hence 

E (J^)'^-^ = E ^*(P)'<^" + E E A.(p)c°„,p</)p„, + E E <n'p5°„"n'</'Pn" • 
i i i n'^p i n'^p 

[n"]^[n'] 

Since n' -< p = [n] if and only if Pn' -< n, and [n"] = [Pn"] ^ [n'] -< n implies that Pn" -< n, we 
can write the above equality in the form 

E (^^°)''^n = E ^^(M)'-^" + E ^"'n'^n' ' (^6) 

This shows that the operator ^,-(J,-°)^ is indeed triangular in the basis B of 7?.„, with eigenvalues 
^j Ai([n])^. It follows from Eqs. (21 ) and (22 1 that H' is also upper triangular in the basis S, and 
that its eigenvalues are given by 



iJ„ = iV52-4^A,([n])2, neN^, 



The latter expression can be simplified by noting that if [n] = p = (pi, . . . ,PAr) and 

Pk-i > Pk = ■■■ = Pk+s > Pk+s+l 
we have 



and therefore 

Afc+,(p) = pk+j + -{l~b)+a{N + j-k-s)= pk+s-j + -{l-h)+a{N + ]-k~s) 
for j ~Q,. . . ,s. Hence 

^ Afe+,(p)2 = ^ (p,+,_^. + -{i-b)+a{N + j-k-s)) = ^ (p, + - (1 - 6) + a(iV - z) 

which yields 

En^Nh'^-^{2p, + l + 2a{N -i)-bf , P = [n] . (27) 



Finally, since H = P-ff'p ^ (cf. Eq. (21 )), the previous discussion implies that H is upper triangular 
in the basis pB of pR-m and that its eigenvalues are also given by Eq. (27) with n G N^. 



C. Spectrum of W 

As discussed at the beginning of this section, the scalar Hamiltonian H^'^ is equivalent to its 
extension H'^^ to the symmetric space A'"^(L^(IR^)), which coincides with the restriction to the 
latter space of the operator H (see Eq. ([l4|)). Since the eigenfunctions of H span the finite- 
dimensional subspace pTZn C L^(M^), for the purposes of computing the discrete spectrum of W'^ 
we can restrict ourselves to the corresponding subspace A^'^{pTZn)- We can construct a basis of the 
latter space by extracting a linearly independent set from the system of generators A^'^{pB), where 
B is the basis of TZn considered in the previous subsection. In this way we easily obtain the basis 
whose elements are the functions 

V'p^pA^^^p^pA^^dle'^'"')' Pi>P2^---^PN, p, e{0,l,...,n}, 

i 

ordered in such a way that ijjp precedes ^p> whenever p ^ p'. It is straightforward to show that 
the operator H'^'^ is upper triangular in the above basis, with eigenvalues given by 

-Ep =iV62-^(2p, + l + 2a(iV-i)-6)^ ^ Pn ^ ■ ■ ■ ^ Pi ^ n . (28) 

i 

Indeed, taking into account that H commutes with the symmetrizer A*^*^ and acts triangularly on 
the basis pB, we have 

iJ^>p - A^'HipcPp) = A^^ {Ep p^p + J2 Cnp PK) = Epi^p + Y, Cnp A'^^(p0„). 

Since 



/ , i-np HVa) — -tJp f p -r ^/^ 
nXp n^p 



we finally obtain 



A-(p^„)=A-(p</.[„])=V[n] 



W''^p=Epi,p+Y.[ H CnpjVp', 



p'^p ^n;[n]=p' 



as claimed. 
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D. Spectrum of H 

The computation of the spectrum of H proceeds along the same lines. Note, first of all, that 
H is equivalent to its extension H to A(L^(M^) (g) E), which in turn is equal to the restriction of 
i? (g) II to the latter space. As before, in order to compute the spectrum of H we should restrict 
ourselves to the finite-dimensional subspace A{pTZn <8) S). A basis of the latter space is obtained 
by extracting a linearly independent set from the system of generators A{pB (g) E). This is easily 
seen to yield the spin functions 

*p,s-pA((/)p|s)), (29) 

where the quantum numbers p and s satisfy 

i) Pi>P2^ ■■■ ^Pn, Pj € {0,l,...,n}, (30a) 

ii) p,=p,, r<3 =^ hl'^' ' = \ (30b) 



The Hamiltonian H is then upper triangular in the basis consisting of the functions (29 1- (30 1, 
ordered so that ^pg precedes vE'p',s' whenever p ^ p'. Indeed, a calculation similar to the one in 
the previous subsection shows that 

#vi/ps = (if® I)(pA(0p|s))) = A(i?(p0p)|s)) 

= A(£;p p(/)p|s) + Y^ C„p /90„|s)) = -Bp*p,s + Yl Cnp *n,s ■ (31) 

Although a given pair of quantum numbers (n, s) in the last sum of the previous equation need 



not satisfy conditions (30), it is easy to see that there is a permutation P (depending on both n 



and s) such that P(n) = [n] and -P(s) = s' do satisfy (30). Since ^n,s differs from the basis vector 
*[n],s' at most by a sign, and n ^ p implies that [n] -< p, all the terms in the last sum of Eq. (31 ) 
precede ^p,s- This establishes our claim and shows that the eigenvalues of H, and thus of H, are 
the numbers 

^p.s ^ Nb'^ -Y (2^' + ^ + ^"(^ - - ^)^ (32) 



where the quantum numbers (p,s) satisfy conditions (30 1 



IV. THE CHAIN'S PARTITION FUNCTION 

In this section we shall compute the partition function of the Frahm-Inozemtsev chain (l8])-(10) 



by exploiting its relation with the dynamical models (IT]) and (pi). Indeed, we have seen in SectionjTT] 
that if we set 6 = /3a with P > 2{N—1), and take the limit a — >■ oo, the eigenfunctions of H become 
sharply peaked around the coordinates of the minimum ^ = {£,1, ■ ■ ■ ,£,n) of the potential U in the 
set ([3]). Note that the condition /3 > 2(A^ — 1) guarantees that the inequality (20) is fulfilled for 
sufficiently large a, so that the Hamiltonians H and H^'^ of the spin and scalar dynamical models 
possess a non-empty point spectrum. From Eq. ([6]) and the definition ([8]) of the Hamiltonian of 
the FI chain, it follows that the eigenvalues of H are approximately given by 

£;y~^f + 4afj, a>l, (33) 

where Ef'^ and £j are two arbitrary eigenvalues of H^'^ and H, respectively. The latter formula 
cannot be directly used to compute the spectrum of the FI chain, since it is not known a priori 



which eigenvalues of H and H^"^ combine to yield a given eigenvalue 8j of V.. However, Eq. (33) 
immediately yields the exact formula 

Z{T) = hm ^(^ (34) 

expressing the partition function Z of the FI chain in terms of the partition functions Z and Z^^ 
of H and H^'^. We shall next evaluate Z by computing the large a limit of the partition functions 
Z and Z"^ 
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A. Partition function of H"" 



Consider, to begin with, the partition function of the scalar model. Expanding Eq. ( 28 ) in 
powers of a we get 

Ep = Ea+iaJ2P'dP + 2i-2N) + Oil), (35) 

i 

where 

Eo = Np^a^ - a^ ^ (/3 + 2i - 2Af) V 2a ^ (/? + 2i - 2N) 

i i 

= \ aN{a{N - l)(3/3 - 2N + 1) + 3{/3 - N + 1)) 

is a constant independent of p. From now on we shall subtract from both H^'^ and H the constant 
energy Eq. With this proviso, when a is sufficiently large the partition function Z'^'^ is approximately 
given by 

Z^^(4ar)~ J2 Y[qP'(P+^^-^^) . 

Note that, by condition (Is]), the coefficient of pi in the RHS of the previous formula is strictly 
positive. In terms of the new summation indices 

rii^Pi- pi+i , 1 < i < TV {pN+i = 0) , 

N 

we have Pi = ^ rij , so that 

N j N 

Y^P^iP + 2i-2N)=Y, n,(/3 + 2i~2N)=Y, n, £(/3 + 2i - 2N) = ^ jn,(/3 - 27V + j + 1) . 

i j'^i J — 1 *— 1 i — 1 

Hence 

Z^=(4aT)~ J2 l[q'"^^^-^^+^+''> , 

ni,...,niv^O j 

where again the coefficient of rij is strictly positive on account of ([5|. For finite a, it is not easy 
to evaluate in closed form the sum in the RHS of the previous equation due to the restriction 
ni + ■ ■ ■ + riN ^ n. This restriction is in fact another peculiarity of the present (hyperbolic) model. 



not present in the trigonometric case |30j . However, since n — >■ oo as a — > cxj on account of Eq. ( 19 1, 
taking into account that q < 1 we finally have 

^hm Z'^4aT) = J2 U l'"' ^^"''^+^'+'' = H £ «'"' (^-2A^+J+i) 

ni,...,nN^O j j nj=0 

= l[{l-q^(^))-\ (36) 

j 

where the dispersion relation J-{j) is defined by 

T{j)=j{P~2N + 3 + l). (37) 

B. Partition function of H 

By Eq. ( 32 ) , the partition function of the spin dynamical model (fTl) can be written as 



O^Pjv^'"^Pi^^ 
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wher e Ep is given by Eq. ( 28 ) and dp is the number of spin quantum numbers s satisfying condi- 
tion ( 30b I . Writing the quantum number p as 



p^ {l^l,...,iyi,...,I^r,---,'^r), i^l^r < ■■ ■ < 1^1 i^n, 



we easily obtain 



nrf^r e = i, 



dp = d(k) 



4=1 



(38) 



(39) 



n © , 



e = -1 



.4=1 



Using the asymptotic expansion (35) and ignoring (as in the scalar case) the constant energy _Eo, 
we have 



Z(4ar) 



2^ dp g • 



O^Pw^'"^Pi^^ 



Setting 






(40) 



and using Eq. ( 38 ) we have 



K, 



^Pi(/3 + 2i-2iV) =^iyj ^ (/3 + 2i-27V) = ^zyjfcj(/3-2iV + 2ii:j -fcj + 1) 

i j = l i=Kj_i + l j=l 

Introducing, as before, the variables 

rii ^ ly^ - Ui+i , 1 < z ^ r [rir = v^) , 

r 

in terms of which i^j = ^ n^ , we can write 

r 

^Pj(/3 + 2i-2iV) = Y^ n^kj{l3~2N + 2Kj-kj + l) =YniNi, 



with 



i=l 



iV, = ^A:^(/3-2iV + 2ii:j -fcj + 1) = (/3-2iV + l)Js:i-^/c|+2 ^ kjki 
= {13-2N+1)K,+J2k]+ Y. k,ki=F{K,), 

j = l l^/#J^4 



(41) 



where T is defined in Eq. ( 37 1 . We thus have 



Z(4aT)^ Y dik) Y. IT' 



,J^{K, 



kePiv ni,...,n^_i>0,n^^O i 



where T'at is the set of all partitions of N with order taken into account. Since F{Ki) > for all 
i, and n — > oo as a — > oo, we finally obtain 



^lim Zi4aT) = ^ <k) E EH '^"'"^"'^ = E <k) E ^"^"^"^^ H E 

kG^N ni,...,nr-i>Onr-^0 i kSPw nr=0 i=lni=l 

r-1 



E <k)(i 



-„^(^-)ri 



ke-PjN 



Hr 



..^■(if.) 



9 



•^(^.) 



n,J=-(A:i) 



(42) 
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Partition function of the FI chain 



Substituting Eqs. (36) and (42) into the freezing trick relation (34) we obtain the following closed 



formula for the partition function of the FI chain: 



N-r 



-,^(K) 



(43) 



kG-Pft 



i=l 



where Ki is given by Eq. (40), r is the number of components of k, and 

.,K'^^,}^{l,...,N}\{Ki,. 



{K[, 



In fact, since Kr ~ 'Y^ ki — N we have 



■,Kr}. 



{K[,..., K'^_^} = {1, . . . , TV - 1} \ {i^i, . . . , Kr-i] . 

It should be noted that the expression (|43l for the partition function of the FI chain coincides 
with the corresponding ones for the HS |30| and PF 44J chains, provided that the dispersion 
relation (371 is replaced by 



-^Hs(j)=j(^-J), -Fpf(j)=j 



From Eq. ( 43 ) it readily follows that the energies of the FI chain are of the form 

£{8) 



JV-l 

E 

1=1 



5,T{i} , 



(44) 



where 5 = [5i- ■ ■ 5n-i) and 6i £ {0, 1}. In fact, Eq. (43) and the discussion in |3?1 imply that the 
numbers 5i are the well-known motifs introduced in Ref. [lOj. More precisely, from Refs. 43, 44] it 
follows that the motifs (with their respective multiplicities taken into account) of the ferromagnetic 
(e = 1) su(m) chain are generated by setting 



|0, Kj+l S^ Ki, 

where the numbers Ki (i ~ 1,...,N) are independent and take the values 0,1, 
particular, in the case of spin 1/2 (m = 2) Eq. (451 is equivalent to 



. . ,rn 



(45) 



1. In 



Si = Kj+l(l- Kj) 



K, e {0, 1} 



(46) 



which together with Eq. ( 44 ) yields the empiric formula for the spectrum proposed in Ref. [27j . 



We shall finish this section by finding the relation between the spectra of the ferromagnetic and 
antiferromagnetic FI chains. From Eqs. (Is]) and (ITT]) it follows that 

(47) 



with 



H^+H^P 


= 2 ^ h,j , 




i^3 


^sinh-2(^^. 





(48) 



The RHS of Eq. (47) is the maximum energy of the ferromagnetic chain when m ^ N, correspond- 



ing to states completely antisymmetric under spin permutations. This maximum energy can be 
easily computed from Eq. (44) noting that J^{i) > for all i, and that the motif 6 = (1 1 • • • 1) is 
compatible with (45) when m^ N (take, for instance, Ki = i). We thus have 



N-l 



2^/1,, 



J2 -^W ^l^iN- l)(3/3 - 4iV + 2) 



(49) 



i=l 



From Eqs. (44), (47) and ( |49| it also follows that if (5i • • • Sn-i) is a motif for the ferromagnetic 
1 — (5jv-i) is a motif for the antiferromagnetic one, and vice versa. This 

for the original HS chain, is a manifestation of 



chain then (1 — Si 

property, empirically discovered by Haldane 

the boson-fermion duality recently established in Ref, 
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V. THE CHAIN'S SPECTRUM 



Over the last few years, there has been growing evidence of the singular character of spin chains 
of Haldane-Shastry type in connection with a number of well-known tests of integrability versus 
chaos in quantum systems. In this section we shall analyze whether the FI chain ([8]) — which is 
most likely integrable [27] — behaves in this respect as other integrable chains of HS type previously 
studied in the literature. To this end, we shall take advantage of the explicit formula (43 1 for the 



partition function to compute the spectrum for relatively large values of N and fixed m — 2,3, ... . 

To begin with, on account of Eqs. (47) and ^49] ) we can restrict ourselves to studying either the 

ferromagnetic or the antiferromagnetic chain (^. Due to the form of the degeneracy factors ( 39 ) 



appearing in the partition function, it is far more convenient to deal with the latter chain (cor- 
responding to e = —1), as we shall do in the rest of this section. In the first place, it is clear 
from Eqs. (37) and (44) that when /? is sufficiently large the energy levels are split into disjoint 
"clusters" centered around the numbers 



/^E^' 



(50) 



where the components of the (antiferromagnetic) motif S are defined in terms of the independent 
numbers k^ G {0, 1, . . . , vn — 1} as 



6,= 



, Ki+l > Ki 



1, 



^i+1 



$ Kj . 



(51) 



It is straightforward to show that the centers of the clusters (50) are the numbers f3j, where j is 
an integer ranging from ^N'{2N — to — niN') to ^N{N — 1), N' = [7V/?tiJ denoting the integer 
part of N/m. Indeed, the minimum value of j corresponds to the motif 



(0---0 10---0 ••• io---o; 



= (A^ - 1) mod 



(52) 



while the maximum value is obtained from the motif (1 1 • • • 1). 

In order to avoid the splitting of the spectrum just described, we shall henceforth assume that /3 
is 0{N) (cf. Eq. ^). With this proviso, we have computed the spectrum of the antiferromagnetic 
FI chain for several values of N, m and f3 using the formula (43) for the partition function. The 



first statistical property that we have analyzed (which is essential for the process of "unfolding" 
the spectrum, as described below) is the cumulative level density (normalized to unity) 



Fi£) 



-N 






(53) 



where di is the degeneracy of the energy Si. It is apparent from our results that, for even moderately 
large values of N, F{£) is very well approximated by the cumulative Gaussian law 



G{£)^ 



1 

/27rcr 



-d£' 



1 



rf 



E — n 

72 cr 



(54) 



with parameters /^ and a respectively equal to the mean and standard deviation of the energy (see 
Fig. fl] left). These parameters can be computed in closed form essentially by taking traces of 
appropriate powers of the Hamiltonian ([s]) , with the result (see Appendix Ib] for the details) : 



M 



Vii'^k)"^" 



l)(3/3-4iV-|-2) 



360 V 



1 



N{N - 1) [16N^ - N^{25f3 - 6) 

+ N{10P^ - 35;9 + 26) + (5/3 - 6)(5/3 + 4)] . 



(55a) 



(55b) 



We have also observed that for /? = 0{N) the level density itself tends to a Gaussian distribution 
when the number of particles tends to infinity, as is the case for other spin chains of HS type whose 
spectrum is equispaced [231110]; see Fig. [11 right. 
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F(£),G{£) 





2000 2200 2400 2600 



1400 1600 1800 2000 2200 2400 2600 



FIG. 1: Left: plot of the cumulative Gaussian distribution (541 (continuous red line) vs. the cumulative level 



density ( |53| l (blue dots) for the antiferromagnetic FI chain with A'' = 20, m = 2, /3 — 40. Right: histogram 
of the density of energy levels (normalized to unity) of the latter chain vs. the Gaussian distribution 

_(£-m)^ 

g{£) = r^ e 2<T^ (continuous red line). 



We shall consider next the distribution of the spacings between consecutive levels in the "un- 
folded" spectrum. Recall [33] that in order to compare different spectra in a consistent way, it is 
necessary to apply to the energy levels fmin = i?i < • • • < £l+\ = ^max the unfolding mapping 
£i '-^ Tji = 'q{£i), where ry(f ) is a continuous approximation to the cumulative level density F{£). 
Indeed, it can be shown that the resulting "unfolded" spectrum {r|i\^^l is uniformly distributed 
regardless of the initial level density. In our case, by the above discussion we can take ri{£) as the 
cumulative Gaussian density (54), i.e., 



1 



Vi = 



1 + erf 



' £i - fJ. 
By convenience, we shall consider the normalized spacings 



(56) 



Hvl+1 - Vi) ' 



i = 1,...,L, 



so that {si}f^i has unit mean. The distribution of the spacings Si is a widely used indicator of 
the integrable vs. chaotic character of a quantum system. More precisely, the celebrated conjec- 
ture of Berry and Tabor posits that for a "typical" quantum integrable system the probability 
density p{s) of the normalized spacings Si should be given by Poisson's law p{s) = e~*. By 
contrast, according to the Bohigas-Giannoni-Schmidt conjecture [35], the spacings distribution 
of a quantum system whose classical counterpart is chaotic should instead follow Wigner's law 
p{s) = (tts/2) exp(— 7rs^/4), as is the case for the Gaussian orthogonal ensemble in random matrix 
theory [33]. Interestingly, the spacings distribution p{s) of all the (integrable) spin chains of HS 
type studied so far [231 EHl EU [38H40] follows neither Poisson's nor Wigner's law. More precisely, 
for these chains the cumulative spacings distribution 

P[s) = f p{s')ds' 
Jo 

can be estimated with remarkable accuracy by the formula 



Pis) c 1 



log(- 



where the parameter 



. f 



Vina 



(57) 



(58) 



is approximately equal to the maximum spacing. In fact, in Refs. |31U40| we have shown that the 
approximation (57)- (58 1 is valid^ for any spectrum satisfying the following conditions: 



^ More precisely |40l . Eq. ||57[l holds in the range so ^ •^ ^ *n 



where sq = Sn 



< Sn 
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(i) The energies are equispaced, i.e., £i+i — Ei = 5£ for i = \, . . . ,L. 

(ii) The cumulative level density (normalized to unity) is approximately given by the Gaussian 
law (ImI). 



(iii) fmax ~ M 1 M ~ '^min ^ CT; wherc jjL and a are the mean and standard deviation of the 
spectrum. 

(iv) fmin and £max are approximately symmetric with respect to /i, namely |iS^min+^max~2/i| ^ 
F — F ■ 

*^max *^min- 

In our case, we have already mentioned that when /3 = 0{N) the FI chain satisfies condition (ii) 
above. It is also easy to check that the third and fourth condition are also satisfied. Indeed, from 
Eqs. ([S]), (|48| and (|49]) it immediately follows that 



6 



(59) 



On the other hand, the minimum energy can be computed from Eqs. (37 1 and (44) using the 
motif (52), with the result 



N' 
fmin = -^ W{N' + 1){2N' + 1) - 'Sm{N' + l)(/3 + 1) + 6N{f3 - A^ + 1)1 . 
6 

Since TV' — N/m + 0(1) and we are assuming that f3 — 0{N), the previous equation yields 

fmin=T^(3/3-4iV) + 0(iV2). 

bm 



From Eqs. (55a), (59) and (61) it immediately follows that 



£n 



H, fi- £n 



4 f 1 - —)n^(3P - 47V) + O(N^) 
12 V rn/ 



(60) 



(61) 



(62) 



and thus fmax — ^min — 0{N^), SO that condition (iv) is clearly satisfied. As to condition (iii), it 
suffices to note that a = 0{N^/'^) by Eq. ( |55b[ ). 

Let us now examine the first condition listed above. The analysis is complicated by the fact 
that the spectrum of the FI chain depends on an essential parameter /3 > 2 (TV — 1), in contrast 
to all the chains of HS type whose spacing distribution has been studied so far. According to our 
numerical computations, there are clearly two different regimes. Indeed, if /3 is an integer or a 
rational number with a "small" denominator, the vast majority of the differences 5£i = £i+i — £i 
are equal to a single value 5£ which depends on /3 (e.g., 5£ = \ when /3 is an odd integer, and 
5£ = 2 when /3 is an even integer). Thus, in this case the spectrum of the FI chain is approximately 
equispaced when N is sufficiently large. Moreover, our computations indicate that in this regime 
the differences 5£i ^ S£ are concentrated in the tails of the level density (see Fig. [2]). As shown in 
Ref. [39], these two properties together with conditions (ii)-(iv) above are sufficient to guarantee 
the validity of the approximation (57). We have verified that this is indeed the case for a wide 



range of values of N, in and /3; see e.g. Fig. [3] for the case N — 24, m — 2 and /3 — 50. 



SSi 



5£i 



4 
3. 
. 2 ■ ■ ■ 

t -3 -2 -1 I 12 3 



FIG. 2; Differences 5Si versus d = {£i — pi) jo for the antiferromagnetic FI chain with N = 24, m = 2 and 
/3 = 50 (left) and AT = 18, m = 3 and /3 = 36.5 (right). 
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FIG. 3: Cumulative spacings distribution and its analytic approximation (57 1 (continuous red line) for the 
antiferromagnetic FI chain with A'' — 24, m — 2 and /3 — 50. For convenience, we have also represented 
Poisson's (green, long dashes) and Wigner's (green, short dashes) cumulative distributions. 



Let us next examine the alternative regime in which the parameter /3 = 0{N) is either irrational 
or a rational with a non-small denominator (greater than 3, say). The key difference with the 
previous case is that the (raw) spectrum is not even approximately equispaced, so that condition 
(i) does not hold. This is essentially due to the fact that the energies £i are of the form f3j + k, 
with j,k integers (cf. Eqs. (37l and (44)). As a consequence, among the differences 6£i there is no 
single dominant value, but rather a set of several most frequent values (cf. Fig. |4J left). 

se. 











., 


0.6 
0.5 
0.4 
0.3 
0.2 






. 



























































O.I 


















FIG. 4: Left: differences 5£i versus e^ = {£i — fij/a for the antiferromagnetic FI chain with A'^ — 24, m = 2 
and /3 = 47 + TT (for convenience we have only shown differences S£i < 0.6, which account for over 98% of 
the total). Right: Cumulative spacings distribution for the latter chain, compared with Poisson's (green, 
long dashes) and Wigner's (green, short dashes) cumulative distributions. 



Our numerical calculations suggest that the presence of several dominant energy differences gives 
rise to discontinuities in the spacings distribution p(s), evidenced by the appearance of "cusps" in 
the plot of P{s) (see Fig. I4J right). The derivation of an analytic approximation to P{s) akin to 
Eq. (571 will appear in a forthcoming paper. 



VI. CONCLUDING REMARKS 



In this paper we have studied the hyperbolic su(7tt-) spin CS model of Aj\j^i type introduced by 
Inozemtsev in Ref. [2S] and its associated spin chain of HS type known as the Frahm-Inozemtsev 
chain [271 . We have computed the spectrum of the former spin model for an arbitrary number of 
particles N and internal degrees of freedom m. Using this result and Polychronakos's freezing trick, 
we have derived a simple closed-form expression for the partition function of the FI chain. With 
the help of this expression, we have shown that the energy levels can be written in terms of the 
standard motifs introduced by Haldane et al. for the HS chain [10] , albeit with a different dispersion 
relation. We have performed a statistical analysis of the chain's spectrum, showing in particular 
that the level density is Gaussian when N is sufficiently large and the chain's parameter f3 is 0{N). 
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It would be desirable to find a rigorous analytic justification of this numerical result, as has been 
recently done for several spin chains of HS type whose partition function factorizes ^32) . We have 
also analyzed the distribution of spacings between consecutive levels of the unfolded spectrum, of 
interest in the context of quantum chaos. Our computations indicate that the density of spacings 
behaves in two qualitatively different ways depending on the parameter /3. More precisely, if j3 is 
an integer or a rational with a "small" denominator, the cumulative spacings distribution P{s) is 
approximately given by the "square root of a logarithm" law typical of other integrable spin chains 
of HS type pll EH |39l SQl |47] . On the other hand, for other values of /? the function P{s) presents 
several "cusps" , but still is neither of Poisson's nor Wigner's type. It would be of interest in this 
respect to ascertain whether this behavior of P{s) is shared by the trigonometric chain of BCn 
type, whose spectrum also depends on a parameter |22j . 

Appendix A: Existence and uniqueness of the minimum of [/ in C 

In this appendix we shall prove that the scalar potential U in Eq. m has a critical point in the 
configuration space (Is]) if and only if /? > 2{N — 1). Moreover, if this condition is satisfied then 
there is a unique critical point, which is in fact a minimum. Our proof is an adaptation of the 
argument in Ref. [IS] for the trigonometric Sutherland potential. 

We shall start by expressing the logarithm of the ground state ( 18 ) of the scalar Hamiltonian (pi) 
as 

ip = log p — aW + Wq , 

where Wq = X)i ^i ^'^d 

i i i<j 

Since, by construction, H^'^c'^ = Eoe"^, where 

Eo^N 

the potential V of the scalar Hamiltonian (pi) can be expressed as 

V = (V(p)2 + Aip + Eo=[ViaW + Wq)] ^ + aAW + Ea . (Al) 

Taking into account that 

V - aAW = a^U + 2(3a ^ e-^"^' , 

i 

equating the coefficients of a^ in each side of Eq. ( |A1[ ) we easily obtain 

U = {WWf + Uo , (A2) 



^ a^{N - l)(3/3 -2N + 1)+ 2a{P + l-N)-l 



with 



[/o = ^A^(7V-l)(3/3-2iV + l) 



Proceeding as in Ref. [3H], we next show that the Hessian of W is negative definite everywhere. 
Indeed, let h = (/ii, . . . , Hn) G M^. Since 

d^W ( _2 V"^ —2 \ —2 

-w—R — = -^ik ( 2/3e ""' + 2^ sinh (x^ - Xj)\ + (1 - ?>ik) sinh {xi - Xfe) , 



we have 






2/3^e-./.^ + i E — ^-^ = -2^5:e--/.^ 



2 fV'sinh (xj-Xj) 



i/i 
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which is clearly negative definite in h for all x e C. Since, by Eq. (A2), 



dU^^yJPVi^dW^ (A3) 

dxi ^ dxidxk dxk ' 

it follows that ^ e C is a critical point of U if and only if it is a critical point of W, i.e., if and 
only if 

Y^ coth(^, - ^j) = /3 - TV + 1 - /3e-^«' , i = 1, . . . , iV . 

Setting 

CN-^+l = /?e-2«' , < Ci < • • • < Cjv , (A4) 

we can easily rewrite the previous system as 

r 

From the last of these equations it follows that 

N-l 

/3 = 2(7V - 1) + Cat + ^ -^ > 2{N - 1) , 

so that Eq. (l5]) is a necessary condition for W ^ and hence f/, to have a critical point in C . 

Conversely, if Eq. ([5]) is satisfied then W must have a unique critical point (a maximum) in C . 
Indeed, since 

Vt/ = (TV „ ^ _ 1) ^ .^^ „ ^ ^ e-2-- + ^(x, + a:,) + E log h"'"' - 6"''^' I 



E -^=G-/3 + 2(TV-l), * = l,...,iV. (A5) 



(2iV - /3 - 2) 5] X, - ^ E e-^- + Y: log |e-^- - e'^^^ 



i<j 



it follows that when /3 > 2(iV — 1) the prepotential W tends to — oo both on the boundary of the 
set C and for Xi — >■ ±00. Hence W must have at least a maximum in C. On the other hand, all 
the critical points of VF in C must be maxima, since we have just seen that the Hessian of W is 
negative definite everywhere. Hence W has at most a relative maximum in C, thus establishing 
our claim. 

The previous result implies that if /3 > 2(A^ — 1) then U has a unique critical point ^ in C. To 



ascertain its nature, it suffices to note that from Eq. (A3) we easily have 



oxiOXj ^—' oxigxk oxkOXi 

■' k -^ 

SO that the Hessian of [/ at ^ is the square of that of W. Since the latter Hessian is negative 
definite, this shows that the Hessian of [/ at ^ is positive definite, so that ^ is a minimum. 

Note, finally, that the zeros of the generalized Laguerre polynomial Lj^ are known to satisfy 



Eq. ( A5 1 , cf . Ref . |35] . From this fact and Eq. ( A4 ) it follows that the coordinates £,i of the unique 
minimum of [/ in C are given by Eq. Q. 

Appendix B: Computation of the mean and variance of the chain's energy 

In this appendix we shall compute in closed form the mean /i and the variance a^ of the spectrum 
of the spin chain (18]) in terms of the number of particles N and internal degrees of freedom m. 
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We shall start with the mean energy jj, = m ^ tiT-L. Using the formulas for the traces of the spin 
permutation operators Sij in Ref . [22] and Eq. ( 49 ) we easily obtain 



'^=(l-^)E%-i^(l-^)^(^-l)(3/5-4^ + 2 



(Bl) 



Consider next the variance of the energy 

2 -N : tni2\ 2 



which is independent of e on account of the identity ( 47 1 . Proceeding as in Ref. [30] we readily 
obtain 



a' = 2 1- 



h)j:^ 



(B2) 



iy^j 



The last sum can be evaluated using the procedure of Ref. [31 , as we shall now show. Note, to 
begin with, that 






E 






E 



a 



(C.-Or frAQ-Q) 



Eh 



2 / .-"'J ■ 

i^3 



(B3) 



In order to evaluate the second of these sums, we use the following identity from Ref. 
\2Q 



E 



(0-0)^ 



= -C/ + 2(/3 + 2)C. - (/3 - 27V + 2)(/3 - 2iV + 6) 



Summing over i and using the identities 



(cf. Ref. [3T]) we easily obtain 

E 






TV 



i5^J 



(G-0)2 12 



— (iV-l)(3/3-4Ar + 8). 



(B4) 



(B5) 



Consider next the first sum in the RHS of Eq. (B3). From Theorem 5.1 of Ref. [50], after a long 
but straightforward calculation wc obtain 



i/i 



(C.-C,)^ 



[8iV" - 8iV(/3 + 1) + 2(3/3 + 4)(^ + 2)] ^ Q 

i 

4(/3 + 2) [AN^ ~ AN{p + 1) + /3(/3 + 2) - 18] ^ G 

N{f3 -2N + 2) [SN^ - 4A^2(3/3 + 2) + iV(6/32 + 8/3 - 216) 

- /3(/32 + 2/3 - 108) + 360] . (B6) 



All the sums in the RHS of this equation can be evaluated from Eq. ( |A5[ ). Indeed, multiplying this 
equation by Q and summing over i we have 

E#V = E|4^-E(cf + c,g + c|) 



(2iv-3)^c? + (EcO' = Ecf-(/3-2^^+2)Ec^ 
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and hence, by the identities (B4|, 



J2 Cf - NlN"^ - 2N^{fi + 1) + 7V(3/3 + 1) + PiP'^ - 1)] 



(B7) 



Similarly, multiplying ( |A5[ ) by Q and summing over i we obtain 



E 



^c/-c^ 



■ / ■ ^i si ■ / ■ si v; 



^(cf + cfc + Gc^ + g*; 



«#J 



2(7V-2)^Cf + 2(^C^)(5:C.)=EC-(/5-2]V + 2)EC?. 



so that, by Eqs. (B4| and (B7), 



J2 Ct ^N{N-f3-l) [N'i^p - 2) - N{f3 + l)(3/3 - 2) - /?(/? - !)(/? - 2)] . 



(B8) 



Substituting Eqs. (B7| and (B8| into Eq. (B6| we easily arrive at 



E 



TV 



(C. - or 720 



(iV - 1) [16N^ + N^{6- 25/3) 



+ iV(10^2 _ 35^ „ 454) ^ 25/3^ + 350/3 + 576] . (B9) 



The previous identity, together with Eqs. (B2), (B3|, and (|B5|), finally yield the explicit for- 



mula (55b) for the variance of the energy of the FI chain. 
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